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ABSTRACT: 

IRAS discovered very luminous objects which emit the vast majority of their radiation 
in the infrared wavelength. The energy source of such a tremendous amount of emission is 
not understood. Starburst and active galactic nuclei (AGN) are thought to be the origin of 
their power. X-ray observations are expected to be able to reveal the characteristics of the 
AGN component. However, some are very X-ray quiet, because the AGN is thought to be 
obscured by a large column density. In order to determine the primary luminosity of the 
nucleus, we need an X-ray observation with a wide energy band at least up to several tens 
of keV. We propose to observe NGC 6240, one of the luminous IRAS galaxies, to investigate 
the characteristics of its AGN, which must be obscured by large column density. 


REPORT: 

Improvements in Background Model, Response and Analysis Software for RXTE: 
Previous to 1998, the background model for the RXTE Proportional Counter Array (PCA) 
did not provide an acceptable subtraction for faint sources such as AGN. An intensive study 
of this problem was carried out at the RXTE Guest observer facility from April to October 
1998. The result was the new L7 background model. This model uses a sum of count 
rates from 7 detectors which have adjacent anodes to predict the rate in the top and mid 
layers. Another ingredient is count rate in the HEXTE particle monitor which is used to 
estimate the radiation dosage from the South Atlantic Anomaly (SAA) and modeled with a 
240 minute half-life. Essential to this model is the assumption that no flux from the object 
can be detected in the bottom layer at low energies. Thus this background model can only 
be applied to objects with count rates < 40 counts s _l . This background model performs 
the best when data taken during the 30 minute interval after emergence from SAA passage 
are discarded, and when data taken when electron rate housekeeping parameter electronm 
greater than 0.1 are discarded. 

The response matrix has also been improved recently. Now, response matrices are con- 
structed for each observation. This makes it possible to account for two time-dependent 
effects: the slow drift in gain of approximately 1% over two years, and a slowly increasing 
xenon fraction in the propane layer. Application of this new response improved the spectra, 
especially in the 10-20 keV range where the flux is low and there is a lot of structure in the 
background. 

Another property that has recently been improved is a better estimation of the uncer- 
tainty associated with the background. Previously, the statistics of the background reflected 
the length of the observation. However, the background model is constructed from many 
background observations, the total length of which is very much longer than a single obser- 
vation. Therefore, the uncertainties in the background spectra should be reduced to reflect 
this fact. A script to correct the background spectra was developed. 
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An improvement to the data analysis software has also recently been made with the 
addition of the REX script. This script permits analysis of a large number of observations 
to be done in a nearly automatic, hands-off manor. The use of this script permits analysis 
of AGN monitoring observations to be done quickly and easily. 

NGC 6240: The RXTE data from NGC 6240 were analyzed in conjunction with archival 
ASCA data. The analysis was complicated by the fact that we found evidence for diffuse hard 
emission contaminating the RXTE spectrum. We found that the spectrum at hard X-rays 
could be characterized by a Compton reflection component; a transmitted AGN component 
is also probably present. The paper has been submitted to MNRAS (Ikebe et al. 1999) and 
a copy is attached to this report. 
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ABSTRACT 

The wide-band energy spectrum of NGC 6240 over the range 0.5-200 keV is 
investigated using the RXTE and ASCA data. The RXTE data provide the 
spectrum beyond the ASCA range (0.5-10 keV) with significant detection of 
signals up to 20 keV and the upper limits above 20 keV. The spectrum above 
10 keV is found to be very flat. A strong iron-K emission line discovered in 
the previous ASCA observation is also confirmed with the RXTE PCA. These 
results provide further evidence for the dominance of a reflection component, 
i.e. emission from cool material illuminated by an AGN. By fitting the spectra 
obtained with RXTE and ASCA simultaneously, we satisfactorily modeled 
the AGN spectrum with a Compton reflection component and probably a 
transmitted AGN component penetrating through a thick absorber. The X- 
ray luminosity of the AGN is estimated to be in the range 4 x 10 43 - 6 x 10 44 
ergs/s in the range 2-10 keV, which categorizes NGC 6240 among the most 
luminous Seyfert nuclei. The ratio of the 2-10 keV X-ray luminosity to the 
infrared luminosity, L x ( 2 - 10 keV)/Lm, is 0.01 - 0.1, which implies a quite 
substantial, if not dominant, contribution of AGN to the infrared luminosity. 

Key words: galaxies: individual: NGC 6240 - X-rays: galaxies 
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1 INTRODUCTION 

The IRAS survey (Neugebauer et al. 1984) discovered many ultraluminous infrared galaxies 
(ULIRGs) that emit the bulk of their energy in infrared (IR) photons. Since their bolometric 
luminosity and the number density are as high as those of quasars, ULIRGs are among the 
most energetic objects in the universe. The most fundamental problem yet to be solved is 
the energy source of the extremely intense infrared emission. 

NGC 6240, a gravitationally interacting system with a complex optical morphology (Fos- 
bury k Wall 1979; Fried k Schulz 1983), is a very interesting example of ULIRG. Its bolo- 
metric luminosity reaches 2.4 x 10 12 Lq (Weight, Joseph, k Meikle 1984; z=0.0245 and 
H 0 = 50 km s -1 Mpc -1 are assumed). NGC 6240 is outstanding in several respects. Its H 2 
1 -» 05(1) at 2.121/xm and [Fell] 1.644/rm line luminosities and the ratio of H 2 to bolomet- 
ric luminosities are the largest currently known (e.g., van der Werf et al. 1993). Further, its 
stellar velocity dispersion of 360 km/s is among the highest values ever found in a galaxy 
centre (e.g., Doyon et al. 1994). 

The energy source of the huge IR luminosity is controversial. Many IR spectroscopic 
studies (e.g. Genzel et al. 1998; Ridgway et al. 1994; Rieke et al. 1985; Weight et al. 1984) 
have suggested that main energy source of the IR emission is starburst activity, which is 
presumably a super-starburst induced by a merger of two galaxies (Joseph k Wright 1985; 
Chevalier k Clegg 1985). The ground-based optical spectrum can be classified as LINER, 
and is interpreted as a result of shock heating (Heckman et al. 1987). On the other hand, 
a significant contribution from an active galactic nucleus (AGN) similar to Seyfert galaxies 
was also discovered from IR spectroscopy (DePoy et al. 1986). Another hint of an AGN in 
NGC 6240 is the presence of compact bright radio cores (Carral et al. 1990 but see Colbert 
et al. 1994). HST discovered a core that is excited higher than LINER (Rafanelli et al. 1997). 

X-ray observation provides an important tool for investigating both the starburst and 
AGN activity. The ROSAT (Triimper 1990) observations showed that NGC 6240 is fairly 
bright in the soft X-ray band below 2 keV with a luminosity larger than 5 x 10 42 ergs s _1 
(0. 1-2.0 keV), and the bulk of the X-ray emission is extended in a scale of ~ 25" (Schulz et al. 
1998; Komossa, Schulz, k Greiner 1998; Iwasawa k Comastri 1998). Detailed spectroscopic 
studies with ASCA (Tanaka, Inoue, k Holt 1994) data of NGC 6240 have shown that the soft 
X-ray spectrum can be explained with two thermal components, a cooler component with a 
temperature of 0.2-0. 6 keV and a hotter component of ~ 1 keV with an excess absorption 
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of ~ 10 22 cm -2 (Iwasawa & Comastri 1998). These results show that these soft X-rays are 
most likely originated from thermal processes, which may arise from starburst activities (e.g. 
Heckman et al. 1987). 

On the other hand, in the 3-10 keV band, another very hard continuum with a strong 
iron-K emission feature was observed with ASCA (Mitsuda 1995; Iwasawa & Comastri 1998). 
The observed spectral feature can be accounted for in terms of Compton reflection from 
optically thick material (e.g. Lightman & White 1988; George, Nandra & Fabian 1990), and 
is generally accepted as evidence for the presence of an AGN in NGC 6240. The emission- 
line profile further indicates that a part of the reflector is highly ionized (Mitsuda 1995; 
Iwasawa & Comastri 1998), and both the ROSAT and ASCA spectra were modeled with a 
reflection component from warm material surrounding the AGN (Netzer, Turner, & George 
1998; Komossa et al. 1998). 

As described above, many observational facts of NGC 6240, in particular X-ray spec- 
troscopy, supports the presence of an AGN in NGC 6240, which may account for a significant 
fraction of the huge IR luminosity. However, the intrinsic power of the AGN is still uncertain, 
The hidden AGN could be visible as a strongly absorbed X-ray continuum above 10 keV, 
penetrating through a thick layer of the obscuring matter. In order to detect the direct X-rays 
from the AGN and to estimate its intrinsic X-ray luminosity, we have observed NGC 6240 
with RXTE, which carries the PC A (Jahoda et al. 1996) and HEXTE (Rothschild et al. 
1998) instruments covering 2-250 keV . 

In the spectral analysis, we also utilized the ASCA data covering 0.5-10 keV energy band 
with two Gas Imaging Spectrometers (GIS: Ohashi et al. 1996) and two X-ray CCD cameras 
(SIS: Burke et al. 1991). Simultaneous use of the ASCA data provides more constraints in 
modeling the spectrum. Moreover, because RXTE has no imaging capability, the RXTE 
spectrum may be subject to contamination by nearby sources. In fact, we find from the 
ASCA image that it is the case, and the ASCA data are used to correct the RXTE spectrum. 


2 RXTE DATA 

2.1 Observation and data reduction 

The RXTE observation of NGC 6240 was performed on 1997 Nov 9-11. The data were 
reduced with Ftools 4.1 and 4.2, We discarded the PCA and HEXTE data that were taken 
when the centre of field of view was within 10° of the local horizon, and when it was off 
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NGC 6240 by more than 0.01 degree. We also selected the time period when all the five 
PCUs of PCA were on. The total on-source exposure time after the data screening process 
is 31.4 ksec and 10.3 ksec for the PCA and HEXTE, respectively. All the five PCUs’ top 
layers of PCA are used for the analysis, while one of 8 detectors of HEXTE (#3 in cluster 
B) is excluded because of its failure. We found no significant time variability from the PCA 
background-subtracted light curve in the entire energy band (3-20 keV) as well as in several 
different energy bands. 

2.2 Derivation of PCA and HEXTE spectra 

We estimated the PCA background using the current standard method. * We attempted 
to estimate systematic error associated with the background subtraction by comparing the 
estimated background with the on-source data in the range above 30 keV where the source 
signal in PCA should be negligible due to low detection efficiency. The count rates of both 
agreed with each other within 0.01%. Therefore, we did not introduce any systematic error 
to the PCA background. 

We estimated the HEXTE background using the data taken during the off-target point- 
ings (Rothschild et al. 1998). In order to estimate the HEXTE-background systematics, we 
accumulated spectra from the +1-5° off-pointings and the —1.5° off-pointings, separately, and 
found that the difference between the two background spectra was within ±2%. Therefore, 
we introduced a 2% systematic error to the HEXTE background spectrum. 

Figure la illustrates the PCA and HEXTE energy spectra thus obtained. Significant 
source signals are detected up to ~25keV in the PCA data. However, the PCA data below 
4 keV and above 20 keV were discarded for the spectral analysis below, due to the uncertainty 
of the energy response reported by Gierlinski et al. (1999) and Jahoda et al. (1996). The 
HEXTE spectrum shows significant signals up to 20 keV, while only upper limits are obtained 
above 20 keV. 

2.3 Spectral analysis 

In order to see the overall spectral characteristics, we first fit the spectra with a simple 
power-law model modified by Galactic absorption, i.e. f{E) = KE- r e ~ a where E 

* The pcabackest software was used with the background model files of pca_bkgd_f aint240_e03v02 .mdl and 
pca_bkgd_f aintl7_e03v01 .mdl. 
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Figure X. (a) The background-subtracted spectrum obtained with PCA and HEXTE (crosses), fitted with a single power-law 
model convolved with the energy-response matrices (solid lines), (b) The ratio between the PCA data and the best-fit power-law 
model. 

is X-ray energy, f(E) is photon flux given in units of photons s _1 cm -2 keV -1 , T is photon 
index, K is the normalization factor, <7 ph is the cross-section of photoelectric absorption given 
by Morrison & McCammon (1983), and 7V H ,Gai is the Galactic hydrogen column density set 
equal to 5.8 x 10 20 cm -2 (Dickey & Lockman 1990). Using XSPEC software (ver. 10.0), we 
performed a minimum-chi-square fitting. The best-fit photon-index is 0.51 yielding x 2 /^ = 
210/32 (Figure la). The ratio of the data to the best-fit power-law model (Figure lb) clearly 
shows evidence for an iron K emission line at ~6 keV and a flatter continuum above ~8 keV. 

Then, we modeled the spectrum with two power-law continua (soft and hard) plus a 
line emission, as expressed by f(E) = ( K s E~ Vs + K h E~ Vh + line{Ki,E l ))e~° p h/VHG »'. The 
iron-line central energy (E t ) and the intensity (K t ) are free parameters, and the line width 
is assumed to be zero. The fit result is shown in Figure 2 and the best-fit parameters are 
summarized in Table 1. Notable properties of the best-fit model are very small photon index, 
~0, of the hard pawer-law continuum, and a large equivalent width, 0.79 keV, of the iron 
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Figure 2. The unfolded energy spectrum (detector response was approximately removed) of the PCA and HEXTE data 
(crosses) fitted with two power-law model and an emission line. Each model component is shown with dotted tines and the 
total model is the solid line. 


Table 1 . The RXTE spectral fit 


Parameter 


Value 

Accounts s“ 1 

cm" 2 at 1 keV) 

1.83 x 10“ 2 

r 8 


4 4+ l 4 
1.2 

Accounts s“ 

1 cm -2 at 1 keV) 

4.45 x 10“ 5 

Th 


0 04 +o 09 
U ' U ^-0.10 

Accounts s _1 

cm -2 ) a ^ 

3.86 x 10" 5 

Ei(keV) 


6.34 ±0.09 

X> 


48.0/28 


a ) : The physical line width is assumed to be zero. The equivalent width is calculated to be 0.79 keV. 

line. An absorption edge structure is noticeable at 7-8 keV in the fit residual (Fig. 2), which 
is most likely the iron K-edge. These features are characteristic of the reflected X-rays from 
an optically thick material as has been pointed out by several authors (Iwasawa &; Comastri 
1998; Netzer et al. 1998). X-rays impinging on optically thick matter are photo-absorbed as 
well as Compton-scattered. These processes form a very flat continuum around ~ 10 keV, 
together with the K-absorption edge and K emission line of iron. We fit the RXTE spectrum 
together with ASCA data in § 4 with models including the Compton reflection. 
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3 ASCA DATA 

3.1 Observation and data reduction 

ASCA observation of NGC 6240 was performed on 27 March, 1994. The GIS provides an 
X-ray image with a circular field of view of ~ 50' diameter. The SIS was operated with 
2-CCD mode, covering a much smaller, rectangular field of view of ~ 11' x 22'. Iwasawa & 
Comastri (1998) presented the result of the SIS data only. Here, we analyse both the SIS 
and GIS data. In particular, the GIS data are useful in order to look for any contamination 
sources in the field around NGC 6240. 

We accepted the GIS and SIS data that were taken when the X-Ray Telescope (XRT: 
Serlemitsos et al. 1995) axis was more than 5° above the local horizon, and when the geomag- 
netic cutoff rigidity was larger than 6 GeV /c in order to ensure a low and stable background. 
Additional screening condition that the elevation angle from the sunlit earth is greater than 
25° and 20° was applied to the SISO and SISl data, respectively. As reported by Turner et 
al. (1997), the ASCA data of NGC 6240 show no significant time variation. 

3.2 GIS image 

The GIS image in the full energy band 0.5-10 keV shows a point-like source that coincides 
with NGC 6240 within ~ 1 arcmin systematic error of the satellite attitude. In addition to 
the emission from NGC 6240, the soft Galactic diffuse emission covers the entire GIS field of 
view. The position of NGC 6240, (/, 6)=(20.73, 27.29), is in the middle of an excess-emission 
structure of Loop I, which is found in the ROSAT All Sky Survey (Snowden et al. 1997). 
The energy spectrum of the excess emission is very soft (kT ~ 0.3 keV), and it is significant 
only in the soft energy band below ~ 2 keV. 

On the other hand, at higher energies, we find an extended structure as clearly seen in 
the image for the 2.4-10 keV band shown in Figure 3, in which the background consisting 
of non-X-ray background (NXB) and cosmic X-ray background (CXB) has been subtracted. 
These are estimated from the data obtained in the blank sky observations and data taken 
when the XRT is pointing at the dark (night) earth (for detail see Ikebe et al. 1995; Ishisaki et 
al. 1997). The background was subtracted to derive the image. In Figure 3, notable emission 
is located north-east of NGC 6240, and extends towards north as well as east. Since the 
extension of the emission is much larger than that of the soft X-ray emission of NGC 6240 
detected with ROSAT, we consider it to be an unassociated contamination source. Since the 
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RA(2000) 


Figure 3. The 2.4-10 keV X-ray intensity contour derived with the GIS data is shown in logarithmic scale. After the background 
(NXB + CXB) is subtracted, the image was smoothed with a Gaussian of a ~ 1.25'. The outer most solid circle shows the 
edge of the GIS field of view-. The position of NGC 6240 is marked with *. 

field of view of RXTE/PCA and HEXTE are both ~1 degree FWHM, the emission from the 
contamination source will contribute to the RXTE spectra significantly. Below we examine 
the ASCA spectrum of NGC 6240 and the properties of this contamination source. 

3.3 Spectrum of NGC 6240 

We construct the GIS energy spectrum of NGC 6240 from a circular region of 3 arcmin 
radius centered on the X-ray peak. The spectra from the two GIS sensors, GIS-2 and GIS-3, 
are summed. The NXB + CXB background is subtracted, as has been done for the GIS 
image. The soft excess emission is not subtracted, but it is negligible above 2.0 keV. 

The SIS spectrum of NGC 6240 is extracted within a circle of 3 arcmin radius from the 
SIS-0 chip-1 and the SIS-1 chip-3. The data taken from the two chips are summed together. 
For the SIS background, we accumulate photons from a region where no contamination 
source is present on the same CCD chips with which NGC 6240 was observed. Therefore, 
the soft Galactic diffuse component is also subtracted as a part of the background. 
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The GIS and SIS spectra thus obtained are illustrated in Figure 4. Following Iwasawa k 
Comastri (1998), we fit the ASCA spectrum of NGC 6240 with a two-component model that 
includes a thermal component and an AGN component. The thermal component consists 
of emission from two optically-thin thermal plasmas of different temperatures, where the 
higher temperature component has an excess absorption. The AGN component consists of 
an absorbed power-law continuum and line emission. The model can be written as; 

f{E) = {thml(T c , Z, K c ) + e ~ a ^ thml{T h , Z, K h ) 

-f line(E t , K t ) + e“ r '> hA ' H - AGN K AGt ,E- r } 

x e _ ' JphNH ' Gil , (!) 

where the parameters in parentheses are kept free. For the thermal plasma emission code, 
we employ the Mewe-Kaastra model (Mewe, Gronenschild, k van den Oord 1985; Mewe, 
Lemen, k van den Oord 1986; Kaastra 1992) modified by Liedahl, Osterheld, k Goldstein 
(1995), which is implemented in XSPEC as the MEKAL model. The metallicities, Z, of the 
two thermal components are assumed to be the same. The abundance ratios among different 
elements are fixed to be the fiducial solar values given by Anders k Grevesse (1989). 

With this model, we fit the GIS and SIS spectra simultaneously. In this simultaneous 
fit, the GIS data below 2.0 keV are excluded because of a slight uncertainty in the energy 
scale of the GIS arising from the complex xenon M-edge structure in this range. A good fit 
is obtained as shown in Figure 4 and Table 2, and the results are essentially the same as 
obtained by Iwasawa k Comastri (1998) using the SIS data only. As noted by Iwasawa k 
Comastri (1998), the striking spectral features are a very flat continuum above ~4 keV with 
a photon index of ~ 0 and a strong iron K-line with an equivalent width of 1.2 keV. 

3.4 Contamination source spectrum 

Although the origin of the extended emission near NGC 6240 is unknown, we only need 
the energy spectrum of this contamination source for the purpose of the present work. We 
construct the spectrum from the whole field of view of the GIS excluding the 3-arcmin radius 
circle centered on the X-ray peak of NGC 6240, from which the estimated background (NXB 
+ CXB) is subtracted. In addition, because of an extended outskirts of the XRT point-spread 
function, the spectrum still contains a significant contribution from NGC 6240. Using a ray- 
tracing simulation software, the contribution from NGC 6240 outside the 3-arcmin radius 
circle is estimated and subtracted. Figure 5 shows the GIS spectrum thus obtained. The 
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Figure 4. (top) The background-subtracted SIS spectra (crosses) fitted with sum of thermal component (dotted lines), an 
emission line (dot-dashed line), and a power-law model (dashed line). All the model components are convolved with the energy 
response matrices. The solid line shows the sum of all the model components, (middle) The same as the top panel, but for the 
GIS data, (bottom) The fit residuals for the SIS data (crosses) and the GIS data (crosses with an open circle). 


Table 2. ASCA spectral fit a > 


Parameter 

Value 

K c (cm 5 ) 6 ) 

3.8 x 10 64 

T c (keV) 

0 Q 1 +0.1 1 

u *' 51 - 0.05 

K h (cm 5 ) 6 ) 

3.5 x 10 65 

T h (keV) 

1 Q 7+ 0 * 24 
l -0.10 

Z (solar) 

1 41+11*1® 

N\\ (10 22 cm -2 ) 

0 Q2 4 " 9 ' 14 
U ’ y ^-0.16 

Accounts s” 1 cm' 2 )^ 

2.9 x 10' 5 

Ei(keV) 

6.46+““ 

Nh.aGN (10 22 cm~ 2 ) 

O.otl 4 

Kagn (counts s' 1 cm' 2 keV -1 @lkeV) 

2.0 x 10“ 5 

r 

-0 in 4-0 ' 29 

o iu_ 0 43 

x 2 / 1 ' 

59.4/55 


a) : Throughout this paper, for fitting the SIS spectrum only, an additional absorption with equivalent hydrogen column 
density of 3.0 x iO 20 cm' 2 was artificially applied, in order to account for a systematic overestimate of the absorption column 
density (Dotani et al. 1996; Cappi et al. 1998). 

^ : Emission measure. 

c) : The physical line width is assumed to be zero. The equivalent width with respect to the power law continuum is 1.24 
keV. 
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Figure 5. (top) The unfolded energy spectrum of the contamination source obtained with GIS. The spectrum is fitted with 
sum of a thermal and a power-law component (dotted lines). Solid line shows the total model, (bottom) The fit residuals. 

resultant GIS spectrum contains the total photons from the contamination source within 
the GIS field of view and the Galactic soft diffuse emission. 

As shown in Figure 5, thus obtained GIS spectrum if fitted satisfactorily (y 2 / v - 18.4/25) 
with the sum of a thin-thermal model (MEKAL model) and a power-law model. For the 
thermal component, we obtain the best-fit temperature of 0.23T0.04 keV and the surface 
brightness of l.OxlO -14 ergs s _l cm -2 arcmin -2 (0.5-2 keV), assuming the element abun- 
dances to be solar. We identify this thermal component to be the Galactic diffuse emission, 
since these values are consistent with the ROSAT results (Snowden et al. 1997). Thus, the 
spectrum of the contamination source can be expressed by a power law with the best-fit 
photon index of 1.56±0.23. Its flux is 2.1 xlO -12 ergs s _1 cm -2 (2 - 10 keV), which is ~ 90% 
as high as that of NGC 6240. Therefore, the contribution from the contamination source to 
the RXTE spectrum is quite substantial, and will be taken into account for the analysis of 
the RXTE spectrum in the next section. We note that the contamination source may ex- 
tend further outside the GIS field of view (see Fig. 3), hence the contribution to the RXTE 
spectrum could be higher than the above obtained flux value. 
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4 RXTE + ASCA JOINT FIT 

In this section, we model the 0.5-200 keV energy spectrum of NGC 6240 obtained with 
RXTE and ASCA, taking into account the contamination source detected in the GIS field 
of view. Below 3 keV, the X-ray emission is dominated by the thermal component that 
presumably originates from the starburst activity in the galaxy, and is well described with 
the two temperature model (Iwasawa & Comastri 1998) as shown in § 3.3. The X-ray emission 
from the AGN in NGC 6240 dominates above 4 keV and we perform detailed modeling below. 

4.1 Power law model 

For modeling the energy spectrum of the AGN, we begin with the simplest model, power-law 
+ line, given by equation 1 in § 3.3, The ASCA GIS and SIS spectra of NGC 6240, the PCA 
spectrum, and the HEXTE spectrum were fitted simultaneously. Since the X-ray intensity 
of the AGN may have varied between the ASCA and RXTE observations performed at 
different periods, the normalization factor of the power-law component, Aagn, is left free 
for the ASCA and RXTE spectra, respectively. Since the iron K-line is presumably generated 
by reprocessing of the X-rays from the AGN on an optically thick matter, the line intensity 
will be proportional to the AGN luminosity on the time average. Therefore, we assumed that 
the ratio of the flux of the iron line to that of the power-law continuum, K t /K AGN , was the 
same in the two separate observations with ASCA and RXTE, respectively. On the other 
hand, the normalization factor of the thermal components, K c and K h , are assumed to be 
common to the RXTE and ASCA spectra, since the thermal component which is extended 
should not vary with time. All other parameters are also tied between the two observations. 
In order to take into account the hard contribution source, a power-law model with the 
photon index 1.56 as derived in the previous section is included in the PCA and HEXTE 
spectra with a free normalization factor. We assume that the power law spectrum extends 
beyond 10 keV with the same photon index. 

As shown in Figure 6 and Table 4.1, the fit is not acceptable. As expected, the obtained 
parameters of the thermal component are essentially the same as those given in Section 3.3, 
since the RXTE spectrum has little influence below 4 keV. However, the best-fit power-law 
continuum of AGN is significantly flatter than that given by Iwasawa & Comastri (1998). 
Furthermore, the residuals require an even harder continuum in 8-15 keV band, and also 
suggests a spectral steepening above ~ 30 keV. These features together with a large equiv- 
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Figure 6. The result of the simultaneous fit to the ASCA/GIS, ASCA/SIS, RXTE/PCA, and RXTE/HEXTE with sum of ther 
mal component, emission line, and a power-law model, (top) The unfolded spectral data (crosses) of ASCA/SIS, RXTE/PCA 
and RXTE/HEXTE. Individual fitting model components are shown in dotted lines, while the total is given with the solid line. 
The GIS data are omitted for the display, (middle) The fit residuals for the SIS data (crosses) and for the GIS data (crosses 
with open circle), (bottom) The fit residuals for the PCA and HEXTE data. 


Table 3, Power law model 


Parameter 

Value 

Ki(ASCA)( counts s _l cm“ 2 ) a ) 

2.8 x 1(T 5 

Ei(keV) 

6.44 

Wh.agn (10 22 cm~ 2 ) 

1.1 

r 

-0.15 

Xagn(ASCA)'') 

1.9 x 10" 5 

K AO v(RXTE)/K Al gu(ASCA) 

1.16 

X 2 /" 

131.5/87 


<0 : The physical line width is assumed to be zero. The equivalent width with respect to the power-law continuum is 1.13 keV. 
: In the unit of counts s” 1 cm -2 keV” 1 @lkeV. 

alent width of Fe K-line are consistent with a Compton reflection spectrum. We therefore 
incorporate Compton reflection in the next subsection. 


4.2 Reflection model 

We employ the Compton reflection model developed by Magdziarz k Zdziarski (1995) im- 
plemented in XSPEC as the PEXRAV model. Using the PEXRAV model, we can calculate 
the reflected X-ray spectrum when X-rays of a power-law spectrum illuminate an optically 
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thick layer of material that is predominately neutral except hydrogen and helium. Here, we 
assume that the direct X-rays from the AGN are totally blocked by a thick absorber on the 
line of sight. The parameters that describe the model are: T, photon index of the incident 
power law spectrum with an exponential cutoff at E c \ R = £l/2ir t where Q is the solid angle 
subtended by the optically thick material; n = cosO , where 6 is the angle between the line 
of sight and the normal vector of the optically thick layer; Ap e , the iron abundance of the 
optically thick material. Element abundances besides iron are assumed to be the solar values. 

In the following fitting, E c is fixed at 200 keV (practically a single power law in the 
observed range), is fixed at 0.45 so that the model spectrum is closest to that averaged 
over all viewing angles as described in Magdziarz & Zdziarski (1995). Combining with the 
thermal component and the line emission, the overall fitting model is given by 

f(E) = {thml(T c , Z, K c ) + thml(T h , Z, K h ) 

+ line(E h Ki) +refl{r,I< AGN )}e^ hNH ^ , (2) 

where free parameters are shown in parentheses. The normalization factors of the reflection 
component, K A G n, for the ASCA and RXTE spectra are left free from each other for possible 
time variation of the AGN luminosity. The flux ratio between the Fe K-line and the AGN 
continuum Ki/R A q^s is assumed to be constant. All other parameters are also assumed to 
be the same between the two observations. Another power-law with a photon index of 1.56 
that represents the contamination source is added to the model for the PCA and HEXTE 
spectrum. 

As shown in Figure 7 and Table 4, this model can account for the 0.5-200 keV wide- 
band spectrum satisfactorily. The flux of the contamination source in the RXTE spectrum 
determined from the fit is slightly larger (~ 10%) than that obtained with the GIS. This is 
probably due to an extension of the source beyond the GIS field of view (see Fig. 3). The two 
normalization factors of the AGN component obtained from the ASCA and RXTE spectra 
respectively happened to be very close (A^ XT ^/A" J ^^ i ' = 0.92— 1.26), as were the Fe K-line 
intensities nearly equal in these two observations. The equivalent width of the iron line with 
respect to the reflection continuum is 0.8 keV, which is a typical value for the fluorescence 
line associated with the Compton-reflection. 

However, the best-fit photon index T of the incident AGN spectrum is 1.26T0.13, which is 
unusually small compared to the typical values for AGNs (e.g. Mushotzky, Done, & Pounds 
1993). In this fit, the viewing angle fi = cosd was fixed at 0.45. Even if we allow n to 
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Figure 7. Unfolded spectrum of the RXTE data (crosses) fitted with sum of the thermal component, reflection component, 
and a power-law for the contamination source. Each model component is shown with dotted lines, while the total is illustrated 
with the solid line. Although the ASCA/GIS and SIS data are not shown here, they are also fitted jointly. 


Table 4. reflection model 


Parameter 

Value (fit l) a > 

Value (fit 2)°> 

r 

1 26 +013 

0.98 + “ 5® 

fj, = cosO 

0.45 (fix) 

6)0 08 tao3 

Fe abundance (solar) 

0 47+ 0 - 22 
u - q ' -0.19 

0 72"^ 0-38 
u - 1 z -0.32 

Kagn(ASCA) c > 

2.7x 10 -3 

4.6x 10 -3 

K’l(ASCA)(counts s -1 cm' 2 )^ 

2.5xl0“ s 

2.4x 10 -5 

E,(ASCA)( keV) 


6.44±0.04 

K AG n(RXTE)/K A gn(ASCA) 

loslg 

1UO -0.15 

X 2 /" 

97.2/87 

96.4/86 


a ) : n is fixed or left free in the fit 1 and 2, respectively. 

b ) : The lower and upper limits correspond to the limitation of the model code. 

c ) : In the unit of counts s~ l cm“ 2 keV -1 @lkeV. 

4) : The physical line width is assumed to be zero. Equivalent width with respect to the reflection continuum is 0.88 keV and 
0.84 keV for the fit 1 and 2, respectively. 

vary (Table 4), photon index larger than 1.54 is not acceptable. The intrinsic luminosity of 
the AGN is estimated to be larger than 4.3 x 10 43 ergs/s, the value corresponding to the 
maximum solid angle, i.e. R = Cl/ 2n, of 1.0. 
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4.3 Reflection + Transmission model 

Although, the above reflection-only model gives a satisfactory fit to the RXTE + ASCA 
data, the small photon index derived in § 4.2 still remains to be problem. As shown below, 
this is resolved by adding an absorbed AGN continuum that is transmitted through a thick 
absorber on the line of sight. (Here we do not consider the extent of the absorber, hence we 
assume no scattering of X-rays into the line of sight.) The transmitted AGN component is 
given by 

e -w H ,AGNKh+ffTh)^ AGN £!-r e -B/Bc ^ (3) 

where cr-rh is the Thomson scattering cross-section, and Nh,agn is the hydrogen-column 
density along line of sight. The attenuation by the Thomson scattering becomes important 
above ~ 10 keV as compared to photoelectric absorption. 

Adding the transmitted AGN component given by eq. 3 to eq. 2, we construct the 
reflected- and transmitted- AGN -I- thermal model as; 

f(E ) = {thml(T c , Z, K c ) + e _<TphNH thml(T h> Z, K h ) 

+ line(Ki,Ei)+refl(r,R,E c ,K^ G N) 

+ e -NH,AGn(° P h+°-n>) K AGli E- r e~ E/Ec } 

x g-^ph^H.Gal (4) 

The photon-index T, the cutoff energy, E c (fixed to 200 keV), and the normalization, A AG n, 
are the same for both the reflection and transmitted components. The relative intensities 
between the transmitted component and the reflection component is determined by the solid 
angle, R , which is left free. The metallicity of the elements lighter than Fe and the viewing 
angle /z = cosd are fixed at 1.0 solar and 0.45, respectively. As in the previous fits, including 
the contamination-source component in the RXTE model, we fit the four spectra jointly 
with this model. 

The reflected- and transmitted- AGN + thermal model gives a good fit (Figure 8, Table 5). 
The derived photon index is V- 1.59 toll- which is consistent with the canonical value of 
1.9 (Pounds et al. 1990). The best-fit solid angle of the reflector, R = Q./2n, is 0.51, and 
the absorption column density, Ah,agNi is found to be 1.7 x 10 21 cm 2 . Since R was poorly 
constrained, the intrinsic luminosity of the AGN has a large error range. If we assume that R 
does not exceed 1.0, L x is obtained to be l.ltj;f x 10 44 ergs/s in the range 2-10 keV. When 
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Figure 8. The fit result of the reftection-and transmitted-AGN model. Crosses are the unfolded spectral data. Each model 
component is shown with dotted lines, while the total is given with solid line. Although the ASCA/GIS and SIS data were also 
fitted jointly, only the RXTE data are shown here. 

we leave the viewing angle /i free, the fitting parameters as well as the intrinsic luminosity 
remain essentially the same, if R < 1.0 (Table 4.3). 

So far, we assumed that the matter that obscures the AGN is located only on the line of 
sight. If the absorbing matter covers a significant solid angle viewed from the AGN such as in 
the case of a torus, some of the incoming photons will be scattered into our line of sight. Then 
the model given by eq. 3 would overestimate the AGN luminosity. Matt, Pompilio, & Franca 
(1999) performed a Monte Carlo calculation of the X-ray transmission through spherically 
distributed matter for various column densities. According to their result, the true AGN 
luminosity would be smaller by a factor of ~ 2, if the absorbing matter is distributed in a 
spherical geometry. Consequently, for the same flux from the cold reflector, the true value of 
R should be larger by a factor of 2 than that in Table 5. The column density of the absorbing 
matter is also subject to a slight overestimation. 

The case that a heavy absorber exists only on the line of sight is rather unlikely. On 
the other hand, since the fitting result shows a relatively small absorption for the reflection 
component, the absorbing matter is probably not covering the entire sphere. There two cases, 
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Table 5. Reflection and transmitted AGN model 


Parameter 

Value (fit 1) 

Value (fit 2)°> 

r 

i cq+0.42 
1,0 - 0.26 

i 71 +0.31 
1 ^-0.37 

rl — cosO 

0.45 (fix) 

&)f) 91 +0.04 
Uyi -0.86 

Fe abundance (solar) 


0 77+ 1 - 23 

' -0.47 

k- A gn(asca)‘> 

9.1X10- 3 

1.3x 10 -2 

Wh.agn (x 10 24 cm -2 ) 

1 7+ 1,0 
1 -0.5 

1 7 +0 - 7 

R 

0 M+ 0 - 99 
U,OA -0.41 

n *ia 4*0.06 0) 

0.34— 02I 

/("((counts s -1 cm -2 ) d ) 

2.4x 10“ 5 * 

2.3xl0~ 5 

E,(keV) 

6 44+0.04 
1,UU -0.I5 

6 44^ 0 04 
°- 44 -0.05 

K agn (RXT£)/Kagn(ASCA) 

I - U ' 5 -0.16 

x 2 A 

86.2/85 

84.8/84 


°) : In the fit, n is left free. Maximum value of R is assumed to be 1. 

k) : The lower and upper limits correspond to the limitation of the model code. 

c ) : In the unit of counts s' -1 cm“ 2 keV -1 @lkeV. 

<*) : The physical line width is assumed to be zero. The equivalent width with respect to the reflection continuum is 0.96 keV 
and 0.95 keV for the fit 1 and 2, respectively. 


a cloud on the line of sight are spherically distributed matter, are considered to represent 
two extremes. 

In conclusion, the reflected- and transmitted - AGN model gives the following AGN pa- 
rameters: r = 1.33 — 2.02, Lx(2-10keV) = 5 x 10 43 - 6 x 10 44 ergs/s, R = 0.1—1, and A^h,agn 
= 1.0 - 2.7xl0 24 cm -2 , where the range of parameters represents not only the statistical 
errors but also the uncertainties due to unknown inclination angle of the reflector, and the 
limits for the two extreme geometries of the absorbing/scattering material. 


5 DISCUSSION 

The X-ray spectrum of NGC 6240 is satisfactorily explained with a model consisting of a 
thermal component and an AGN component. For the AGN component, and equally good 
fit is obtained either by the reflection-only model or by the reflected- and transmitted - AGN 
model. Hence, we cannot conclude the detection of the direct AGN component from the 
present data. However, without the direct (transmitted) component the derived photon 
index is noticeably smaller than the typical values for AGN. We consider it more probable 
that the direct AGN component is present at high energies. 

According to the unified model for Seyfert galaxies (e.g. Antonucci 1993), one can inter- 
pret the result in terms of a tilted molecular torus in which the near-side of the torus acts as 
an absorber and the far-side acts as a reflector. If the torus is axisymmetric, both the reflect- 
ing part and the absorbing part would have a similar column density. With the reflected- 
and transmitted - AGN model, the column density of the obscuring matter is estimated to 
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be ~ 2 x 10 24 cm -2 , corresponding to Thomson optical depth of ~ 1. This is practically 
enough for a reflector. Thus, the present result is consistent with the interpretation that 
NGC 6240 has a similar geometry to Seyfert 2 galaxies. The obtained column density gives 
an IR extinction of Ak ~ 100, which is consistent with the fact that the previous IR obser- 
vations did not find the AGN itself. Narrow line region that has not been seen in the optical 
observations would be obscured as well. 

The intrinsic X-ray luminosity of NGC 6240 is estimated in § 4.3 to be in the range 
5 x 10 43 — 6 x 10 44 ergs/s in the range 2-10 keV. For the reflection-only case (§ 4.2), the 
luminosity is estimated to be larger than 4 x 10 43 ergs/s, while the upper bound of luminosity 
is not determined. However, it would be plausible that the solid angle factor R is not much 
less than several %, which gives essentially the same luminosity range as the above. This 
luminosity is among those of the most luminous Seyfert nuclei, and even comparable to those 
of quasars. The bolometric luminosity of the AGN may well exceed 10 45 erg/s. 

Concerning the power source of the huge IR luminosity, a measure of the contribution of 
an AGN to the IR emission is a ratio of X-ray luminosity to IR luminosity, L x /L m . Based 
on our results given above and the IR flux calculated from the IRAS result with the formula, 
Fir = flux(25/xm)x (z^m ) + flux(60pm) x (i/ 60/im ), gives L x (2 — 10 keV)/L [R = 0.01 - 0.1 
for NGC 6240. This agrees with those of other Seyfert nuclei given by Ward et al. (1988). 

The BeppoSAX results on NGC 6240 have recently been published by Vignati et al. 
(1999). The authors claim that the direct AGN component was positively detected, though 
they assumed the photon index of the AGN to be 1.8. Other than that, their results and 
interpretation are essentially in agreement with ours. 

The luminosity of the AGN can be used to infer the mass of the central black hole in 
NGC 6240. If the AGN accretes about 1% of the Eddington luminosity, the total luminosity 
of 10 45 erg/s would lead to a mass of M A gn — 10 9 A/ o , which is consistent with the value 
expected from the galaxy-mass to black-hoie-mass relation (Lauer et al. 1997), with the 
estimated galaxy mass of 10 u-l2 M o after having completed its merging epoch (Shier & 
Fisher 1997). 

6 CONCLUSION 

The 0.5-200 keV wide band energy spectrum of NGC 6240 obtained with RXTE and ASCA 
is accounted for in terms of a soft thermal component and a hard AGN component. The 
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soft component is presumably due to star burst activity. The AGN component consists of a 
Compton reflection component accompanied by an intense Fe-K emission line and probably 
a transmitted component (a direct component penetrating through a thick absorber). The 
detection of the transmitted component is not conclusive from the fitting. However, without 
a transmitted component, the photon index is unusually small. Assuming that the solid 
angle factor R does not exceed 1, we estimated the intrinsic X-ray luminosity of the AGN 
in the range 2-10 keV to be in the range 4 x 10 43 - 6 x 10 44 ergs/s, which yields the ratio 
of the X-ray luminosity (2-10 keV) to the IR luminosity of 0.01 — 0.1. The column density 
that obscures the central AGN is estimated to be larger than 1.0 x 10 24 cm -2 . These results 
show that NGC 6240 is among the most luminous Seyfert 2 galaxies. 
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